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Abstract. It is known that the Pak-Stanley labeling of the Shi hy- 
perplane arrangement provides a bijection between the regions of the 
arrangement and parking functions. For any graph G, we define the 
G-semiorder arrangement and show that the Pak-Stanley labeling of its 
regions produces all G-parking functions. 

In his study of Kazhdan-Lusztig cells of the affine Weyl group of type 
An-i, [15], J.-Y. Shi introduced the arrangement of hyperplanes in M" now 
known as the Shi arrangement: 

Xi — Xj = 0,1 1 < i < j < n. 

Among other things, he proved that the number of regions in the complement 
of this set of hyperplanes is (n + Cayley's formula for the number of 
trees on n + 1 labeled vertices. The first bijective proof of this fact is due to 
Pak and Stanley, [16], who provide a method for labeling the regions with 
parking functions of size n. Given a graph G, Postnikov and Shapiro, [14], 
introduced the notion of a G-parking function. In the case where there exists 
a vertex q connected by edges to every other vertex of G, Duval, Klivans, and 
Martin, [7], have defined the G-Shi arrangement, and conjecture that when 
its regions are labeled by the method of Pak and Stanley, the resulting labels 
are exactly the G-parking functions with respect to q. In this case, however, 
there may be duplicates among the labels. Letting G be the complete graph 
on n -I- 1 vertices recaptures the original result of Pak and Stanley. Our work 
was motivated by this conjecture. 

Figure 1 serves as a guide to our paper. The four corners of the square 
in Figure 1 are labeled by structures we associate to a graph G, which we 
now describe. 

Parking functions and superstable configurations. Chip-hring is a 
key tool for us. In the abelian sandpile model for G, one first chooses a 

vertex of G to serve as the "sink." Then placing grains of sand (or chips) 
on each of the nonsink vertices defines a configuration, c, on G. A vertex v 
is unstable in c if it has at least as many grains of sand as its degree. In that 
case, one may fire (or topple) v by sending one grain of sand from v to each 
of its neighbors. In this process neighbors of v may become unstable, them- 
selves. Sand is not necessarily conserved under firing: grains of sand that 
are sent to the sink vertex disappear. For undirected, connected graphs — to 
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Figure 1. Schematic diagram of results. 



which we limit ourselves here — by repeatedly firing vertices, c is eventually 
transformed into a stable configuration, i.e., one with no unstable vertices. 

Instead of firing one vertex at a time, one may consider a firing rule that 
allows sets of vertices to be fired simultaneously. This gives rise to a stronger 
type of stability and a corresponding set of superstable configurations (see 
Section 2.4). These superstable configurations serve to define G-parking 
functions. A G-parking function with respect to the chosen sink is a function 
from the vertices of G to the integers whose value at the sink is —1 and 
whose values at the nonsink vertices are the numbers of grains of sand for 
some superstable configuration on G. Thus, there is only a slight difference 
between superstable configurations and G-parking functions, and at any 
rate, they are in one-to-one correspondence. Parking functions are more 
widely-known, and they originally arose independently (for instance in [11, 
p. 545]). This accounts for why they, and not superstable configurations, 
are mentioned in the title of this paper. However, given the importance of 
chip-firing to our work, we will mainly refer to superstable configurations 
outside of this introduction. 

Although G-parking functions and superstable configurations arc defined 
with respect to a chosen sink vertex, we have found it natural and convenient 
to first prove our results in a sinkless context, akin to working in projective 
space rather than affine space. To this end, we define generalized G-parking 
functions, which we call quasi- superstable divisors on G. These appear in 
the bottom left corner of Figure 1, denoted by S. As indicated by the figure, 
quasi-superstables are formally defined in Definition 10. 
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Semiorders and the semiorder arrangement. The semiorder arrange- 
ment, [17], is the set of n{n — 1) hyperplanes in M" given by 

Xi - Xj = 1, i,j G {1, . . . ,n}, j. 

Its regions are in bijection with certain n-element posets called semiorders. 

A G-semiorder is a semiorder whose elements are the vertices of G. The 
collection of G-semiorders, denoted by I, appears in the bottom right corner 
of Figure 1. In the same way that a Shi arrangement is modified by Duval, 
Klivans, and Martin in [7] to take into account the structure of a graph, we 
modify a semiorder arrangement to produce the G-semiorder arrangement 
associated with G. The regions of the G-semiorder arrangement are denoted 
by TZ in the top left corner of Figure 1. 

Fixing a sink vertex, v, we then refine the definition of a G-semiorder 
arrangement to get the (G, v) -semiorder arrangement (Definition 28) . 

Orientations. A partial orientation of G consists of orienting some, not 
necessarily all, of the edges of G. The G -semiorientations, denoted by O 
in the top right of Figure 1, are the partial orientations satisfying an extra 
condition relative to the cycles of the underlying graph, G. 

Relations among the structures. Using the method of Pak and Stanley, 

we label each region of the (G, t;)-semiorder arrangement with what amounts 
to a function from the vertices of G to the integers taking the value — 1 
at V. The main goal of this paper is Theorem 32, which, along with 
Theorem 34, shows that the set of Pak-Stanley labels that are negative only 
at the sink are exactly the set of G-parking functions. The set of Pak-Stanley 
labels consists of only G-parking functions exactly when the sink is adjacent 
to every nonsink vertex. 

Our main goal, just described, is a consequence of the analogous result in 
the nonsink context, represented by the left- most vertical arrow in Figure 1. 
This arrow represents the fact that the Pak-Stanley labels for the regions 
of the G-semiorder arrangement form the set of quasi-superstables, which is 
proved in Theorem 23. This theorem is proved by establishing the relations 
represented by the other arrows in Figure 1. 

In working on the G-Shi conjecture of Duval, Klivans, and Martin, we 
were led to labeling the regions of the G-Shi arrangement with partial ori- 
entations of G. We used Farkas' lemma to develop a rule determining which 
partial orientations would appear as labels. However, it seemed that the cri- 
terion we developed was awkward, and that it would become much simpler 
and symmetric if we altered the hyperplane arrangement somewhat. This 
is how we were led to consider G-semiorder arrangements. The resulting 
correspondence between the regions of the G-semiorder arrangement and 
their labels with partial orientations (the G-semiorientations) is represented 
by the top horizontal arrow in Figure 1 and is established in Theorem 14. 

The correspondence between maximal superstable configurations — those 
containing the most amount of sand — and acyclic orientations of the edges 
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of G has been noted several times, in different forms ([3], [4], [8], [9]). In [3], 
the correspondence is a consequence of an extended version of Dhar's burn- 
ing algorithm, a tool from the chip-firing literature, [6]. The input to their 
extended algorithm is a maximal superstable configuration, c. The config- 
uration is transformed into an unstable configuration c. The sequence of 
vertex firings that stabilizes c is then used to orient the edges of G. We 
modify this algorithm in Section 3.3 to apply to all superstable configura- 
tions, not just the maximal ones. This algorithm is encoded (in its nonsink 
version) as the bottom horizontal and right vertical arrows in Figure 1 . 

Ultimately, we label each region of the G-semiorder arrangement with (i) 
a G-semiorientation, (ii) semiorders on the vertices of G, and (iii) a quasi- 
superstable divisor on G. Our aim — to show that all quasi-superstables 
appear as labels — is achieved by showing that the upper-left and the lower- 
right triangles in Figure 1 commute and that the mapping referred to as 
"indeg — 1" along the diagonal in Figure 1 is surjective, (cf. Theorem 23, 
Corollary 20, and Theorem 19, respectively). 

Organization. Following this introduction, the paper is organized into four 
main sections and a conclusion. Section 1 is an extended example that il- 
lustrates our main results and may serve as a foundation for understanding 
them. The rest of the paper may be considered a justification of the claims 
made there. Section 2 defines our four main graphical structures. Sec- 
tion 3 contains the central results describing the correspondences between 
these structures, as discussed above. Section 4 explains how to transfer the 
results of Section 3 to the context in which a sink vertex is chosen. The con- 
cluding section presents a conjecture and suggests further lines of inquiry. 
Imagine starting with the G-semiorder arrangement but then perpendic- 
ularly displacing the hyperplanes (replacing the original hyperplanes with 
parallel ones) before applying the labeling method of Pak and Stanley. As 
one slides the hyperplanes, some regions disappear and new regions form. 
We conjecture that as long as a "central region" is preserved, the resulting 
set of labels does not change. A special case implies the G-Shi conjecture of 
Duval, Klivans, and Martin. 

Acknowledgments. We thank Art Duval, Caroline Klivans, and Jeremy 
Martin for encouraging us to work on the G-Shi conjecture and for helpful 
comments. We thank Collin Perkinson for help with proofreading. We also 
especially thank the anonymous referee, who helped to improve the quality 
of our exposition. 

1. Introductory example 

We introduce our main results with an example, beginning with an ex- 
planation of the construction of Figure 3. Consider the arrangement of six 
planes in M^, 

Xi-Xj = l, i,j e {1,2,3}, i ^ j. 
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The complement of these planes in consists of 19 connected components. 
The planes are all parallel to the vector (1, 1, 1), so for our purposes it suffices 
to intersect the arrangement with the perpendicular plane, xi + ,T2 + .T3 = 0, 
as pictured in Figure 3. Our goal here is to explain the labeling of the 19 
regions by vertex-labeled, partially oriented graphs. 

To start, each region contains a copy of the graph, G, pictured in Figure 2 
after removal of the sink vertex, q. The vertices of the remaining triangle are 



V3 




Q 



Figure 2. Graph G. 

labeled by integers, and some of the edges of the triangle are oriented. We 
refer to the vertex labels as vectors, (01,02,03), where Cj is the label for Vi, 
and designate oriented edges as ordered pairs of vertices, {u,v), where u is 
the tail and v is the head. If there were more room, each region would be 
labeled by the full graph, G. In each region, the label on q would be — 1, 
and each edge incident on q would be oriented with tail at q. 

The labeling of the regions is done inductively, starting at the center. 
In our example, the central region of the arrangement is a hexagon, and 
its triangle has vertex labels (0,0,-1) and no oriented edges (besides the 
assumed ones from the sink). The rule for the central region is that a 
vertex Vi has a label of if {q,Vi} is an edge; otherwise the label is — 1. 
More generally, the rule for any region is that the vertex Vi is labeled by one 
less then the number of oriented edges pointing into Vi. 

The central region shares an edge with six bordering regions. Let r denote 
one of these regions, and say that Xj — = 1 is the border between it and 
the central region. Moving from the central region to r, one may think of 
the value of xj as increasing at the expense of Xi. We record this fact by 
orienting the edge {vi.Vj} as {vi,Vj) and increasing the vertex label for Vj 
by one. Label the five remaining bordering regions similarly. For example, 
moving into the compact region directly above the central region, we cross 
into the region where > X2 + 1- Thus, the label for this compact region 
adds an oriented edge {v2,vs), and the label for ^3 increases from —1 to 0. 

At this point in the labeling process, there would be six unlabeled regions 
bordering the six that were just labeled. Let r be one of these unlabeled 
regions, and let — = 1 be the edge it shares with a region r' that was 
just labeled. In our case, there are two choices for r' , but this choice does 
not affect the eventual label for r. To label r, start with the label for r', add 
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X-\_ — X2 "hi ^3 — Xi -|- 1 Xi — Xo ~h 1 X2 — X^ -|- 1 




Figure 3. Labeled (G, g)-semiorder arrangement. 

the oriented edge (xfc,x^), and increase the vertex label for vi by 1. After 
labeling these six regions, there are six remaining unlabel regions, and these 
are labeled by continuing the procedure just described. 

In the end, each region is labeled with a copy of G having labeled vertices 
and a partial orientation of its edges. The resulting partial orientations have 
a special property. Call an unoriented edge "blank." Consider a collection 
of edges, C, forming a cycle in G. Under any partial orientation, some 
of the edges in C will have an orientation. If it is possible to orient the 
remaining blank edges in C to get a directed cycle, then the special property 
is that there must be a greater number of blank edges in C than oriented 
edges: "more blanks than arrows for potential directed cycles." For example, 
consider the label for the region just above the central region. It has one 
oriented edge, (^2,^3). By orienting the two blank edges as (^1,^2) and 
{vs,vi), we would get a directed cycle, but this potential cycle has two 
blanks and only one arrow. Theorem 14 guarantees that the 19 regions are 
in bijection with partial orientations of G having this special property. 
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As mentioned above, the vertex labels are given as one less than the 
indegree at each vertex. In our example, there are eight distinct vertex 
labels with nonnegative entries: 

(0, 0, 0), (1, 0, 0), (0, 1, 0), (0, 0, 1), (2, 0, 0), (1, 0, 1), (0, 2, 0), (0, 1, 1). 

The surjectivity of the map ip in Theorem 19 ultimately implies, through 
Theorem 32, that appending —1 for the sink, q, to each of these labels yields 
the set of G-parking functions (with respect to q). As they are, these labels 
are exactly the superstable configurations for the abelian sandpile model 
on G. Subtracting each from the maximal stable configuration, (2,2,1), 
gives the recurrent configurations, i.e., the elements of the sandpile group: 

(2, 2, 1), (1, 2, 1), (2, 1, 1), (2, 2, 0), (0, 2, 1), (1, 2, 0), (2, 0, 1), (2, 1, 0). 

Now pick any region r and a point t = {ti,t2,t^) G r. The point deter- 
mines a collection of closed intervals, li = [ti,ti + 1] . Turn the set of intervals 
into a poset, P, by saying li < Ij if Ij lies completely to the left of Ij with 
no overlap, i.e., if tj + 1 < tj. Overlapping intervals are not comparable 
in P. Posets arising from finite sets of intervals, ordered in this way, are 
called semiorders. Identifying interval Ij with vertex Vi gives a poset on the 
nonsink vertices, and we then set q to be the unique minimal element to get 
a poset on all of the vertices. 

For instance, the point p = (1, 3, 2.5) is in the unbounded region directly 
to the right of the central region in Figure 3, giving the collection of intervals: 



The Hasse diagram for the resulting semiorder on the vertices of G is shown 
in Figure 4. 




Figure 4. Semiorder determined by (1,3,2.5) G M^. 

The partial orientation of the label for region r can then be read from the 
semiorder: {vi,Vj) appears as an oriented edge if and only if {vi,Vj} is an 
edge of G and Vi < Vj in P. In general, varying the point p selected in r 
may result in different semiorders: they may disagree for pairs of vertices 
that do not determine an edge of G. (In our example, each pair of the Vi 
form an edge, so only one poset on the vertices arises from each region.) 
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How does this example generalize to an arbitrary graph, G? Suppose G 
has designated sink q and has nonsink vertices vi,...,Vn- Form the arrange- 
ment of hyperplancs, xi — Xj = 1 for all i ^ j such that {vi,Vj} is an edge 
of G. Thus, for each edge, {vi,Vj}, there is a "stripe" consisting of the two 
hyperplanes Xi — Xj = ±1. The partial orientations in our labeling record 
on which side of each stripe a given region lies. 

Label the central region, for which |xj— Xj| < 1 for all edges {vi, Vj}, by the 
graph G, orienting the edges incident with q so that they point away from q. 
Label the Vi by the number of oriented edges pointing into Vi minus 1 (hence, 
by or —1), and label q with —1. Then proceed inductively to label all of 
the regions, as in the example. In the end, the regions will be in bijection 
with those partial orientations satisfying the property that for each potential 
directed cycle, there are more blanks than arrows. The collection of vertex 
labels with nonnegative values at each Vi are the G-parking functions, which 
correspond to the elements of the sandpile group. Picking a point in any 
region determines a semiorder on the vertices of G, from which one may 
reconstruct the labeling of G for the region. 

We now describe a version of the above construction that avoids an initial 
choice of a sink. Let G be a graph with vertices {vq, . . . ,w„}. Form the 
hypcrplane arrangement as described above: Xj — Xj = 1 for each i j such 
that {vi, Vj} is an edge of G. Thus, there are two hyperplanes for each edge. 
Label the central region — given by |xi — Xj| < 1 for all i,j such that {i,j} is 
an edge of G — with a copy of G having no oriented edges and with a —1 at 
each vertex. Proceed as before, labeling each region. By Theorem 26, those 
regions for which Vi has label —1, and all other vertices have nonnegative 
labels, are the G-parking functions with respect to Vi. 

For example, again take G to be the graph in Figure 2 but with no vertex 
chosen as sink (take q = vq). The corresponding hyperplane arrangement 
consists of ten hyperplanes in M*^. Each hyperplane is parallel to the vector 
(1, 1, 1, 1), so we project onto the hyperplane given by Ylf=o — 0; which we 
identify with R^, to get an arrangement whose 109 regions are in bijection 
with those of the original arrangement. The central region is a polytope 
with ten faces. It is pictured in Figure 5 along with two of its bordering 
regions. The bordering region forming a pyramid on top of the central region 
is labeled by a copy of G with one directed edge, {v2,vi), and with vertex 
label (— 1, 0, — 1, — 1). The other bordering region in the figure is labeled 
by a copy of G with one directed edge, (u3,ui), and with the same vertex 
label, (— 1, 0, — 1, — 1). Figure 6 depicts those (unbounded) regions of the 
hyperplane arrangement for G that satisfy Xi > xq + 1 for i = 1, 2, 3. They 
are in bijection with the regions in Figure 3, and corresponding regions 
would have the same labels. (Recall that for convenience the drawing of the 
sink vertex and its edges is suppressed in Figure 3.) 
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Figure 5. The central region and two bordering regions for G. 




Figure 6. Regions for G corresponding to choosing a sink vertex. 

2. Four structures associated with G 

From now on, we take G to be a finite, connected, undirected graph, 
with vertices V = {fo, . . . , Vn} and edges E. Loops and multiple edges are 
disallowed (the former for convenience of notation). 

2.1. G-semiorientations. In this section we define O, the collection of G- 
semiorientations. A partial orientation of G is a choice of directions for a 
subset of the edges of G. Formally, a partial orientation is a subset O C 
V xV with the property that if (u, v) E O, then {u, v} ^ E and (u, u) ^ O. 
Let O be a partial orientation. If e = {m, v} ^ E and (u, v) G O, then 
despite the ambiguity, we write e G O and say e is oriented. In that case, we 
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think of e as an arrow from u to v and write = u and e'^ = v. If neither 
{u,v) nor {v,u) is in O, we write e ^ O and say that e is an unoriented 
or blank edge. The outdegree of the vertex u £ V relative to O, denoted 
outdego(M), is the number of edges e G O such that e~ = u. Simiharly, the 
indegree oi u E V relative to O, denoted 'uidegQ{u), is the number of edges 
e E O such that e+ = u. We use the notation deg{u) to denote the ordinary 
degree of u, i.e., the number of e € E' containing u. Fixing O, some of the 
edges of any cycle C C. E oi G will be oriented and others will be blank. If 
it is possible to assign directions to the blank edges so that C would become 
a directed cycle, then we call C a potential cycle for O. 

Definition 1. A G-semiorientation is a partial orientation, O, such that 
each potential cycle for O has more blank edges than oriented edges. The 
set of G-semiorientations of G is denoted O. 

2.2. G-semiorders. A reference for ordinary semiorders is [17]. Let k be 
any nonnegative integer, and consider a collection of unit length closed in- 
tervals, P = {Ii, . . . , Ik}, of the real line. Order the elements of P by /j < Ij 
if Ii lies strictly to the left of Ij, i.e., if Ij = [oj, + 1] and Ij = [aj,aj + 1], 
then ai + 1 < aj. Any poset isomorphic to a poset P, constructed as above, 
is called a semiorder. The number of non-isomorphic semiorders with k 
elements is the k-th. Catalan number, Ck, and there is a corresponding gen- 
erating function 



If fk denotes the number of labeled semiorders, there is the exponential 
generating function 



Definition 2. A G-semiorder is a semiorder on the vertices of G. The set 
of G-semiorders is denoted I. 

2.3. The G-semiorder arrangement. 

Definition 3. The G-semiorder arrangement, denoted y, is the set of 2\E\ 
hyperplanes in M"'*^-'^ given by 



for all i ^ j such that {vi,Vj} G E. The regions of ^, denoted TZ, are the 
connected components of M""^^ \ J^. 

If G were the complete graph on n + 1 vertices, the G-semiorder arrange- 
ment would be the ordinary semiorder arrangement discussed in [17], whose 
regions are in bijection with labeled semiorders on n + 1 elements. For gen- 
eral G, each region of the G-semiorder arrangement is a union of regions 
from the ordinary semiorder arrangement. 





k>0 
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2.4. Quasi-superstables on G. Now designate q := vq as the sink vertex. 
We recall the basic facts about the abelian sandpile model on G, including 

G-parking functions. A reference for the sandpile results stated here is [10]. 
(More references for the abelian sandpile model: [2] and [5] are seminal; [6] 
and [13] are general references; and [12] is a quick overview.) We then 
proceed to define 5, the set of quasi-superstable divisors of G. 

Let D = diag(deg'yo, . . . , degi;„), and let A be the adjacency matrix for G, 
defined by 

'■^ 1 otherwise. 
The Laplacian matrix for G is 

A = D-A. 

The reduced Laplacian is the matrix A obtained by removing the first row 
and column of A, i.e., the row and column corresponding to the sink vertex. 

A configuration on G is an clement of the free abelian group on the nonsink 
vertices. Having ordered the nonsink vertices, as above, we identify the set 
of configurations with Z" in the natural way: c = Y17=i Ci '"i (ci, • . . , c„). 
A divisor on G is an element of the free abelian group on all of the vertices 
of G, which we similarly identify with Z""*""^. Once a sink is chosen, we may 
consider configurations as those divisors whose sink coefficient is 0. 

Given two configurations or two divisors c and c', we write c > c' if Cj > 
for all i. We say c is nonnegative if c > 0, i.e., if each component of c is 
nonnegativc. 

If X is a subset of the nonsink vertices, write Ix for the configuration 
whose i-th component is 1 if G X and otherwise. Firing X from a 
configuration c results in the configuration c — A Ix- If X = {vi}, we call 
this operation firing Vi. One may speak of firing the sink vertex, which 
adds 1 to each vertex connected to the sink. Since the sum of the columns 
of the Laplacian matrix is zero, firing the set of all nonsink vertices is the 
same as reverse-firing the sink, subtracting 1 from each vertex attached to 
the sink. 

Definition 4. Let c be a nonnegative configuration on G. We say c is 
stable if there is no i such that firing Vi from c results in a nonnegative 
configuration. We say that c is superstable if there is no nonempty subset X 
of the nonsink vertices such that firing X from c results in a nonnegative 
configuration. 

The notion of a superstable configuration is essentially the same as that 
of a parking function. 

Definition 5. A function /: V ^ Z is a G-parking function (with respect 
to q) if f(q) = —1 and (/(f i), . . . , f{vn)) is a superstable configuration on G. 
We identify a G-parking function / with the divisor ^^^y f{v) v. 
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Suppose c is a nonnegative configuration. Then c is stable exactly when 
Ci < degwj for all i. A vertex Vi such that Ci < degwj is said to be stable in c; 
otherwise it is unstable. Firing a set of vertices is legal from c if the resulting 
configuration is nonnegative. In particular, firing a single unstable vertex 
of c is legal. A sequence of vertices is called a legal firing sequence for c if 
each vertex in the sequence is unstable after firing the previous vertices in 
the sequence. Since there is a path from each nonsink vertex to the sink 
in G, there is a legal firing sequence leading to a stable configuration c° 
called the stabilization of c. This process is called stabilizing c. It turns out 
that c° is independent of the order in which unstable vertices are fired, as is 
the number of times each vertex is fired in reaching c°. 

Definition 6. A stable configuration c > is recurrent if given any nonneg- 
ative configuration a, there exists a nonnegative configuration b such that 

c = (a + b)°. The recurrent elements with the operation of (vertex- wise) 
addition followed by stabilization is called the sandpile group of G (with 
respect to q), denoted Sand(G). 

It is well-known that the sandpile group actually is a group and the map- 
ping 

Sand(G) ^ Z"/image(A) 
c !->■ c 

is an isomorphism. Each equivalence class of modulo the image of the 

reduced Laplacian contains a unique recurrent element. It is also known 
that each equivalence class contains a unique superstable element. Define 
the maximal stable configuration to be 

n 

Cmax = '^{degVi -l)vi. 

i=l 

The next two propositions are well-known. 

Proposition 7 ([10, Theorem 4.4]). The configuration c is recurrent if and 
only if Cmax — c is superstable. 

Proposition 8 (Dhar's burning algorithm, [6], [10, Lemma 4.1]). Let b>0 

be a stable configuration on G, and let b be the configuration obtained from b 
by firing the sink. The following are equivalent: 

(1) b is recurrent, 

(2) {b)° = b, i.e., the stabilization of b is b, 

(3) in stabilizing b, each nonsink vertex fires exactly once. 

Remark. We refer to Proposition 8 as Dhar's burning algorithm, although it 
would more properly be called the theoretical underpinning of the algorithm. 

Definition 9. Let K{G) be the graph G with the addition of a new vertex q 
and edges {q, Vi} for all vertices Vi, including q = vq. Set q as the sink vertex 
oiK{G). 



ORIENTATIONS, SEMIORDERS, ARRANGEMENTS, AND PARKING FUNCTIONS 13 



Thus, the set of nonsink vertices of K(G) is V, and each of these is 
connected to the sink, q, by an edge. The divisors on G are exactly the 
configurations on K{G). 

Definition 10. A divisor c G Z"^^ on G is called quasi- super stable if c = 
c — ly for some superstable c on K{G). The collection of quasi-superstable 
divisors is denoted S. 

For the relation between the superstables and quasi-superstables of G, see 
Theorem 26 (2). 

3. CORRESPONDENCES BETWEEN THE STRUCTURES 
So far, we have defined the following structures on G: 
O : G-semioricntations, 
I : G-semiorders, 
^ : the G-semiorder arrangement, 
S : quasi-superstable divisors on G. 

In this section, we describe relations among these structures, culminating in 
Theorems 19 and 23. 

3.1. Semiorders and semiorientations. 

Definition 11. A G-scmiordcr P and a G-semiorientation O are compatible 
if for each edge e = {u, v} of G, we have that u < v \{ and only if {u, v) € O. 
Thus, if e ^ O, then u and v are not comparable in P. 

Given a G-semiorder, P, define 

Op = {{u, v) : {u, v} £ E and u < v in P}. 

Theorem 12. Let P be a G-semiorder. Then Op G O and Op is the unique 
element o/O compatible with P. 

Proof. The only part of this theorem that is not immediate from the defi- 
nitions is the fact that every potential cycle for Op has more blanks edges 
than oriented edges, which we now prove. Let a = {ei, . . . , e^} be a poten- 
tial cycle. We may assume that = {■Ui,«i+i} where itfc+i = ui, and that 
for each i, either (i) (itj,'Ui+i) G Op, in which case < ttj+i, or (ii) is 
a blank edge, in which case Ui and are not comparable. Since P is a 
semiorder, it is isomorphic to a semiorder of intervals, allowing us to identify 
each Ui with an interval Ij = [oj, + 1]. If Cj G Op, we have /j < li+i, in 
which case Oj < Oj+i — 1; and if Cj ^ Op, then Li and J^+i overlap, so in 
particular, Oj < Oj+i + 1. Thus, if there are S oriented edges and /3 blank 
edges in a, since Uk+i = ui, 

ai + 5 — (3 < ai, 

with equality if and only if /3 = (5 = 0. However, since a has at least one 
edge, the inequality must be strict, and 13 > 6 as required. □ 



14 



SAM HOPKINS AND DAVID PERKINSON 



Thus, we can refer to the semiorientation determined by a G-semiorder 
and define the mapping 



V4 




Vl V2 



Figure 7. The house graph. 

Example 13. Let G be the house graph of Figure 7. Figure 8 depicts a 
G-semiorder and its corresponding compatible G-semiorientation. □ 

3.2. A bijection between semiorientations and hyperplane regions. 

We now define a mapping 

If C G O, let p{0) be the region defined by the following inequalities: for 
each edge e of G: 

• if e = {vi, Vj) G O, then Xj > Xi + 1, 

• if e ^ O, then \xi — Xj\ < 1. 

The reader may find it helpful to skim the example below the proof of 
the following theorem before reading the proof itself. 



V4 




P unit intervals Op 

Figure 8. A G-semiorder P on the vertices of the house 
graph, G; a collection of intervals realizing the semiorder; 
and the unique G-semiorientation, i^{P) = Op, compatible 
with P. 
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Theorem 14. The mapping p is a well-defined bijection. 

Proof. Let O G O, and let r = p{0). The system of inequalities defining r 
indicates on which side of each hyperplane of the region sits. To show r 
is a region of it suffices to show that r is nonempty. 

We define a directed, weighted graph G' with the same vertices as G in 
the following manner: for each edge e between vertices of G, 

• if e = (wj, f j) G C, then (uj, Vj) G G' and the weight of (uj, Vj) in G' 
is —1, 

• if e = {vi^Vj} ^ O, then {vi,Vj), {vj,Vi) G G' and the weight of both 

{vi,Vj) and {vj,Vi) in G' is 1. 

Choose an ordering of the edges of G': e'^, 62, . . . , e'^. Define a k x (n + 1) 
edge- vertex adjacency matrix with rows ri, . . . , as follows: if = {vi, Vj), 
let r£ be the vector having 1 in the ith entry, —1 in the jth entry, and Os 
elsewhere. Let 6 be a column vector in M'^ where bi is the weight of e^. 

Thus, the inequalities of p{0) are encoded as Ax < b. By Farkas' lemma 
the insolvability of Ax < 6 is equivalent to the existence of a row vector 
y = {yi,-,yk) satisfying: 

(1) 2/i>0 Vz, y^O, yA = 0, yb<0. 

For sake of contradiction, suppose such a y exists. The support of y is 

supp(y) = {i:yi^ 0}. 

Among all row vectors satisfying condition (1), suppose y has been chosen 
so that the cardinality of its support is minimal. Say ii G supp(y) and 
^'h ~ Hence, r^^ has a —1 in its Vj-th entry. Since 

yA = yiri \-ykrk = (i, 

and the components of y arc nonnegativc, there must be some £2 G supp(|/) 
such that r^j has a 1 in its vj-th entry. This row will have a —1 in some 
other entry, forcing the existence of some £s G supp(y) such that r^g has 
a 1 in that entry, and so on. Since the support of y is finite, the sequence 
^1,^2; • • • , eventually has a repeat. Thus, there is a sequence of elements 
ji := im+i,j2 ■= £m+2, ■ ■ ■ , jt ■= ^m+t in the support of y for some t corre- 
sponding to a directed cycle of edges e'j^,. . . , e'j^ in G'. 

Let z = (zi, . . . , Zk) be the row vector with 2;^ = 1 if £ G {ji, . . . , jt} and 
zi = 0, otherwise. Since the support of z corresponds to a directed cycle of 
edges in G', we have zA = 0. Furthermore, since any potential cycle in p{0) 
has more blank edges than oriented edges, we have z ■ b > 0. Let a = 
minjyjj, . . . ,yjt} and define y' = y — az. Then y' satisfies condition (1) but 
its support is strictly contained in the support of y, yielding a contradiction. 
So there must be some solution to Ax < b, which means that r = p{0) G TZ. 

We now define a mapping 



T : 7^ ^ O. 
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If r G 7?. let r(r) be the partial orientation of G obtained by (i) (vi, Vj) G r(r) 
if {vi, Vj} e E and Xj > Xi + 1 in r, and (ii) all other edges of G are blank. 
Once we show r is well-defined, it is immediate that it is the inverse of p. 

Let r & TZ. To see that T(r) G O, suppose r(r) has a potential cycle a 
having at least as many oriented edges as blank edges. Define A and b 
as above to encode the system of inequalities that defines the region r as 
Ax < b. Let y be the row vector with Is in the entries corresponding to the 
(oriented and blank) edges of a and Os elsewhere. We have y > 0, y 0, 
yA = 0, and y ■ b < 0. But by Farkas' lemma this means Ax < b has no 
solutions, contradicting the fact that r eTZ. Thus, T(r) GO. □ 

Example 15. This example is intended to be read in conjunction with the 
proof of Theorem 14 and illustrates Farkas' lemma at work. Consider the 
partial orientation, O, of a triangle pictured in Figure 9. Naively attemt- 

A 

Figure 9. Partial orientation, O. 



ping to apply the region-labeling function, p, to O results in the system of 
inequalities Ax < b where 



A 



1 

1 

-1 



-1 
1 








-1 

-1 
1 



Xi 
X3 



The rows of A correspond to directed edges, {vi,V2), {v2,V3), (vijV^), and 
(^3, Vl), respectively, of the underlying triangle, G. The last two inequalities 
combine to describe the "sandwich" \xi — X3I < 1, corresponding to the 
undirected edge, {i;i,i;3}. 

The partial orientation, O, has a potential cycle with more oriented edges 
than blanks, and hence is not a G-semiorientation. By the way in which A 
is defined in the proof of Theorem 14, we know that adding the rows of A 
corresponding to the directed cycle {vi,V2), {v2,V3), (f3,fi) gives the zero 
vector. Our system of inequalities is inconsistent since adding the corre- 
sponding entries of b gives —1 (since there are more oriented edges than 
blanks in O for this cycle). □ 

Let r be a region in the image of p. t = {to, . . . ,tn) & r, define the unit 
intervals Ii = [ti,ti + 1] for i = 0, . . . ,n. Define Pt to be the G-semiorder 
determined by these intervals, labeled by the vertices of G by identifying 
with Vj. 
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Pt Pf 

t = (0, 1.5, 2.75, 3.5, 2) t' = (0, 2, 4.75, 5.5, 4) 

Figure 10. Two semiorders. 

Theorem 16. Let O G O, and let r = p(0). The semiorders Pt as t ranges 
over points in r are exactly the G-semiorders compatible with O. 

Proof. Choose any t e r. Say {vi, Vj) G O. Then Xj > Xi + 1 in r; so 
tj > ti + 1, and hence, Vi < Vj in Pf. Now suppose e = {vi,Vj} G E but 
e ^ O. Then \ti — tj\ < 1, which means that li and Ij overlap, and hence, Vi 
and Vj are not comparable in Pt. This shows that Pt is compatible with O. 

Now let P be a G-semiorder compatible with O. The scmiordcr P is 
isomorphic to the semiorder on a set of unit intervals, {Ii}2=o, where li 
corresponds to Vi. Say /j = [ti,ti + 1] for each i, and let t = (to, . . . ,ti). 
So P = Pt, but we must show that t e r. Suppose e = {vi,Vj} G E. If 
Xj > Xj + 1 in r, then since O = r(r), we have {vi,Vj) G O, and thus 
tj > ti + 1. li \xi — Xj\ < 1, then e ^ O and the intervals /j and Ij overlap, 
i.e., \ti — tj\ < 1. Hence, f G r. □ 

Example 17. Let G be the house graph pictured in Figure 7, and let O be 
the G-semiorientation pictured on the right in Figure 8. The region p{0) 
corresponding to O is determined by the system of inequalities: 

X2> Xi + 1, X3 > X2 + 1, 

X4 > X5 + 1, X5 > Xi + 1, 

|X3 — X4I < 1, |X3 — X5I < 1. 

Figure 10 displays points t, t' G p{0) whose corresponding semiorders, Pt 
and Pf, are the two G-semiorders compatible with O. □ 

3.3. Superstables algorithm. So far, we have described relations among 
three of our four structures: G-semiorders, G-semiorientations, and the re- 
gions of the G-semiorder arrangement. We would now like to forge a con- 
nection between these and the quasi-superstables on G (ultimately relating 
all of this, in Section 4, to ordinary superstables). 

Dhar's burning algorithm says that starting from a recurrent configura- 
tion, then firing the sink, each non-sink vertex will fire exactly once in the 
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stabilization process. We use these vertex firings to build a partial orienta- 
tion starting from the unoriented graph. To sketch the idea, suppose that 
at some point in the stabilization process there are nonsink vertices v and w 
such that e = {v, w} is an edge and v is unstable. Suppose e has not al- 
ready been oriented or marked as a blank edge. When v fires, if w is stable, 
orient e as {v,w); otherwise mark e blank. Since all vertices fire, each edge 
is visited. 

Starting with a quasi-superstable, c, on G, Proposition 7 provides a cor- 
responding recurrent configuration on K{G). We use Dhar's algorithm as 
above to create a G-semiorientation, O, of G. Using the construction of O 
as a guide, one may further refine the procedure to simultaneously create 
a G-semiorder compatible with O. This is the idea behind the intervals 
labeled by J in the algorithm, below. Theorem 19 then provides the im- 
portant connection: subtract one from each entry of the indegree sequence 
relative to O to recapture c. 

A careful description of our algorithm follows. The input is a quasi- 
superstable divisor on G and a vertex ordering, and the output is a G- 
semiorder P and the G-semiorientation Op compatible with P. Letting 6 be 
the symmetric group on {0, . . . , n}, identify cr G © with the vertex ordering 
(Vo-(o)! • • • ,Va{n))- Then our algorithm defines two mappings : 5 x 6 — I 
and : 5 X 6 — 7- O. Theorem 19, to follow, validates the algorithm and 
shows that it produces a commutative diagram 

(2) sxe -^-^ I 

o, 

where v is the mapping defined in Section 3.1. 

Given a c G 5 and a vertex ordering cr G 6, the algorithm proceeds as 
follows: 



INITIALIZATION 

Let Cmax be the maximal stable configuration of K{G), and let 

b — Cmax C. 

Since c G 5, we can write c = c — ly for some superstable c on K{G). 
Therefore, b = (cmax ~ c) + ly, i.e., b is the configuration on K{G) obtained 
from Cmax — c by firing the sink q of K{G). By Proposition 7, Cmax — c is a 
recurrent configuration on K{G)\ so by Proposition 8, every element of V 
will fire exactly once while stabilizing h. 

Let ui,...,Uk be the vertices that arc unstable in b. Take this list of 
vertices to be ordered according to a, that is, if Ui = v^(^£-j and uj = v^^rn) 
with £ < m, then i < j. For each Ui, associate an interval, 

J{Ui) = [i/{k + l),l + i/{k + l)]. 
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Thus, all the J(ui) overlap. Form a queue, 

Q = {ui, . . .,Uk). 

Initialize the partial orientation of G as O = 0. 
LOOP 

Repeat the following until the queue is empty: 
Say the queue is Q = {wq, . . . , Wi), and 

J{wo) = [a,a + 1]. 
If £ > 0, define e using the interval for wi: 

J{wi) = [a + £,a + £ + 1]. 

Otherwise, take e = 1. 

Fire wq and replace b by the resulting configuration. Remove Wq from the 
queue and mark it so that it will never again appear in the queue. 

Let zi,. . . ,zt be the vertices that just became unstable with the firing 
of and that have not yet been fired by the algorithm, listed in order 
according to a. Define 

J{zi) = [a + 1 + i£/(t + l),a + 2 + i£/(t + 1)]. 

Add the to Q, in order: 

Q = {wi,...,we,zi,...,zt). 

For each edge e = {wq, v} that is not oriented or marked as blank, 

(1) if V was already unstable before the firing of wq, mark e as blank, 
so e will not subsequently be added to O; 

(2) otherwise orient the edge out from v, i.e., add {wq,v) to O. 

OUTPUT 

The intervals {J{v) : v G V} determine a semiorder. Identifying v 
with J{v) gives a semiorder, P, on V. Define (/)(c) = P and r]{c) = O. 
□ 

Example 18. Let G be the house graph of Figure 7, and fix the vertex 
ordering, wi, . . . , ^5. Then c = (—1, 0, 0, 0, 0) is a quasi-superstable divisor 
on G. Figure 11 illustrates the application of the superstablcs algorithm 
to c and the given vertex ordering, producing a G-semiorientation, O. The 
dotted lines emanating from each vertex to the exterior of the house graph 
represent the edges to the sink of K{G). 

We have Cmax = (2,2,3,2,3) for K{G). The algorithm starts with the 
configuration h = Cmax — c in the top left corner of Figure 11. The vertex vi 
is unstable in b, and firing produces two oriented edges and two unstable 
vertices, V2 and v^. Since V2 comes first in the vertex ordering, it is fired next. 
The subscripts keep track of vertex-firing precedence. Note that when v 5 is 
fired, the edge {^3,^5} is marked blank since, by that time, 'U3 is unstable. 
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Proceeding clockwise around the diagram, the algorithm terminates at the 
bottom left corner. 

Let c := c+ly, a superstable configuration on K{G). When the algorithm 
terminates, the resulting configuration is Cmax — c, which is recurrent by 
Proposition 7. Adding ly reproduces the starting configuration c. However, 
adding ly is equivalent to firing the sink of K(G), which explains (via Dhar's 
burning algorithm) why every vertex was guaranteed to fire. □ 
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Figure 11. Superstables algorithm (cf. Example 18). Start 
in the upper-left corner. Bold numbers connote unstable ver- 
tices. Subscripts indicate firing precedence. 



Theorem 19. Given P G I, for each vertex v let np{v) denote the cardi- 
nality of the set 

{u E V : u < V and {u, v} G E}. 

There are mappings 



^{np{v) - l)v, 



and 



O y^(indeg0(f ) - 1) v. 
vev 
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Let TT : 5 X © — )■ 5 be the first projection mapping, and let (p and r] be the 
mappings defined by the superstables algorithm. 
The following diagram commutes: 

(3) 5x6 




All mappings in the diagram are surjective except possibly for and r]. 

Proof. The image of (p is in I. In stabilizing b, eacti vertex eventually 
becomes unstable and is assigned an interval. Thus, the semiorder produced 
by the algorithm is a semiorder on all the vertices of G. 

The image of 9 is in S. Let Pel. Then P is isomorphic to the semiorder 

on a collection of intervals, {J{v) : v G V}. Letting c = 9{P), we must 
show that c + ly is a superstable configuration of K{G). Let Cmax be the 
maximal stable configuration on K{G), and define b = Cmax — (c + ly). Note 
that 6 > 0. Starting with b and firing the sink of K{G) gives b = Cmax—c- By 
Propositions 7 and 8, we must show that there is an ordering of V forming 
a legal firing sequence for b. 

Choose any a e & so that a{i) < a{j) if min Jj < min Jj. In particular, 
this means that if Vi < Vj in P, then a{i) < a{j). The v-th component of b 
is 

K = dcg^(G.)(f) - np{v). 

We can legally fire the vertices in the order given by a since when it be- 
comes v's turn to fire, it will have received np{v) grains of sand from the 
firings of those vertices u neighboring v such that u < v and will thus be 
unstable. 

We have 6 o (j) = n. Let P = <p{c, a) and b = Cmax — c where Cmax is 
the maximal stable configuration on K(G). Let v & V . If u is unstable 
in 6, then = — 1 and there are no vertices smaller than v in P. Hence, 
np{v) — 1 = Ct, = — 1, as required. Otherwise, in the course of the algorithm, 
say u is the vertex whose firing causes v to become unstable. Then the 
vertices that are less than v in P are exactly the vertices w such that w < u. 
Among these w, only those that are attached to v by an edge contribute to 
making v unstable. Thus, exactly np{v) grains of sand are added to f in 6 to 
make v unstable. Since by = <^QgK{G) v — l — Cy, we have that np{v) = 1 + 
as required. 

The mapping v is surjective. The surjectivity of v follows from Theo- 
rem 16. 
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We have o v = and the image of i/; is in S. Let P £ I and 
O = iy{P). Since O is compatible with P, we have that np{v) = 'm.degQ{v) 
for each v £ V. Hence, ij){v{P)) = 6{P) as mappings of configurations. 
Since v is surjective, the image of i/j is contained in the image of 6, hence 
in S. 

We have i/ o ^ = 77. Given (c, cr) G 5 x (S, let O = 77(0, a) and P = (l){c, a). 
We must show that O is compatible with P. Let b be as in the algorithm. 
Run the algorithm up until it is a vertex u's turn to fire. Suppose {u, v} G E. 
Then v being stable at this point is equivalent to u < v in P and equivalent 
to {u,v)eO. □ 

Corollary 20. Let a e &. Then the mappings 

defined by (paic) = 4>{c,a) and rjaic) = r]{c,a), are injective with left in- 
verses 6 and tjj, respectively. 

Let 5max denote the maximal quasi- super stables under the relation "< " 
defined in Section 2.4, and let Omax denote the acyclic orientations of G 
(elements of O in which each edge is oriented). For each c E S, define the 
degree of c to be deg(c) := ^„gy c„- Then, 

(1) the restriction of ija to iSmax gives a bijection 

(2) for c & S, we have deg(c) < g — 1, where g := \E\ — \V\-{- 1 is the 
genus of G, with equality if and only if c £ 5max ■ 

Proof. Both (pa- and rja are injective on S with the claimed left inverses since 
the restriction of tt in diagram 3 of Theorem 19 to 5 x {a} is bijective. 

To show T?cr(5max) ^ Omax, let C G 5 and suppose that O := r]a{c) ^ Omax- 

According to Theorem 16 there is a point t = {tQ,...,tn) £ M"'"'"^ and a 
corresponding G-semiorder, Pt, such that iy{Pt) = O. For ease of notation, 
we may assume that ^ " " " < in- Choose t' G M"+^ such that t'- > ti for 
all i and t'-^^-^ > t'^ + 1 for < i < n - 1. Let O' := v{Pt,) and c' := ^{O'). 
Then the set of oriented edges O is a proper subset of O' . Hence, c < c', 
showing c ^ <Smax, as desired. 

Now let c G 5, and choose any O G O such that c = ip{0). From the 
definition of tp, we have deg(c) = \0\ — \V\ < \E\ — \V\ = g — 1, with eqiiality 
exactly when O G Omax- If c < c' for some quasi-superstable c', then since 
deg(c) < deg(c'), it follows that O ^ Omax- This shows, for instance, that 

V'(Omax) Q 5max- 

It remains to be shown that ip restricted to Omax is injective. Let c G 5max 
and choose any O such that c = '4>{0). We have seen that O G Omax- From 
the definition oiip, we see that indeg0(t;) is determined by c for each v £V. 
Then outdeg0(f) = degQ^v) — indeg^ since O orients every edge of G, and 
as noted in [3], these outdegrees determine O. To see this, consider the 
graph G oriented by O. Since O is acyclic, this oriented graph must have 
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Figure 12. Figure for Example 22. 

sinks, i.e., vertices with outdegree 0. Thus, c determines the orientation of 
all edges incident on these sinks. Now remove these sinks and their incoming 
edges. The orientation of the resulting graph is still acyclic, and its sinks 
may also be determined from c. Iterate to see that O is determined by c, 
and hence, i\} is injective when restricted to Omax- 

Theorem 21. The mapping v is bijective if and only if G is a complete 
graph. 

Proof. The surjectivity of v is part of Theorem 19. Suppose that G is a 
complete graph and that z^(Pi) = i'{P2) for some P\,P2 € I- Then Vi < Vj 
in Pi implies {vi,Vj) £ v{Pi), which in turn implies Vi < Vj in P2. Similarly, 
Vi < Vj in P2 implies Vi < Vj in Pi. Thus, Pi = P2, so v is injective. 

Now suppose G is not complete, i.e., there exist Vi,Vj G V such that 
{vi,Vj} ^ E. Let Pi be the G-scmiorder where no vertices are comparable, 
and let P2 be the G-semiorder where Vi < vj and all other pairs of vertices 
are incomparable. We have z^(Pi) = = i^{P2)- So is not injective. □ 

Example 22. The mappings (p and rj are not necessarily surjective. Con- 
sider the case where G is a cycle graph on five vertices. Figure 12 displays the 

Hasse diagram of a G-semiorder, P, and its compatible G-semiorientation 
z/(P) = Op. For any ordering a of the vertices, r]{tl;{Op),a) = O' where O' 
is as pictured in Figure 12. □ 

3.4. Labeling regions with superstables. We can inductively label the 
regions of the G-semiorder arrangement, J^, with quasi-superstablcs. We 
call this labeling the (generalized) Pak-Stanley labeling after [16]. Start 
with the central region of the region defined by \xi — Xj\ < 1 for all i,j 
such that {vi,Vj} is an edge of G. Label this region with the divisor — ly. 
Put the central region in a queue, Q. Then, for as long as Q is not empty: 

(1) Remove the first region r from Q. 

(2) For each unlabeled region r' bordering r: 

(a) Determine the unique indices i ^ j such that \xi — Xj\ < 1 in r 
but Xj > Xi + 1 in r' . 

(b) If r is labeled by c = J2k=o ^k^k-, then label r' hy d = c + Vj. 
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(c) Add r' to the end of Q. □ 

For each r e TZ, define A(r) to be the label assigned to r by the above 
algorithm. Theorem 23 guarantees that when the algorithm terminates, 
the regions are labeled with quasi-superstables and the labels are, in fact, 
independent of the order in which regions are removed from the queue. 

Recall the bijection p: <0 ^ TZ, and let r = Also recall the mapping ^ 
from Theorem 19 and the mapping r/^ from Corollary 20. 

Theorem 23. We have X = ip or. So there is a surjective mapping 

A: TZ 

and, for each a & &, a commutative diagram 
(4) 




Proof. All the statements in the theorem follow directly from the equality 
X = o T, which we now prove by induction. If r is the central region 
then r(r) is the partial orientation in which all edges are marked blank. 
Therefore, in this case, A(r) = ilj{T{r)). 

Let the labeling algorithm run, and suppose that the region r has just 
been removed from the queue. By induction, suppose that X = tp o t when 
restricted to those regions that have been labeled so far. Let r' be an un- 
labeled region bordering r. Say that the inequalities that define r are the 
same as those that define r' except that \xi — Xj\ < 1 in r and Xj > Xi + 1 
in r'. It follows that if O = T{r) then T(r') = OU {{vi,Vj)}. Therefore, 
going from T(r) to T(r'), only the indegree of Vj has increased by one, so 

A(r') = A(r) + Vj (algorithm) 
= ip{T{r)) + Vj (induction) 
= V(T(r')). 

The result follows by induction. □ 



4. Fixing a sink vertex 

In this section, we finally fully explain Figure 3. For that figure, we 
started with a graph with a given sink vertex and constructed a hyperplane 
arrangement with regions labeled by partial orientations of the graph and 
sandpile configurations. The nonnegative configurations arising were exactly 
the superstables for the graph. 

Let G be a graph with vertices V = {vq, . . . , Vn}, and designate vertex vq 
as the sink. 
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Definition 24. The {G,vo)-semiorientations, denoted Oq, are the G-semi- 
orientations satisfying the additional requirement that vq is a source: 

Oo = {O G O : outdego(uo) = deg{vo)}. 

The set of admissible {G, i;o)-semiorientations is 

Oo = {O G Oo : indego('yi) > 1 for all i ^ 0}. 

Thus, while Vo is the sink for the sandpile model on G — i.e., the sink 
for the sake of defining the sandpile group, superstable configurations, and 
G-parking functions — it is a source for any O € Oq. 

Example 25. Figure 3 shows the admissible (G, g)-semiorientations for the 
graph in Figure 2. Due to lack of space, the requisite edges {q, vi) and {q, V2) 
are understood in Figure 3, but not drawn. □ 

Our next goal is to describe the image of Oq under the mapping tp: O ^ S 
from Theorem 19, thus accounting for the configurations labeling the regions 
in Figure 3. We show below that the image of Oq is the set of quasi- 
superstables assigning the value —1 to vq and whose only other negative 
values must occur at vertices not connected to vq by an edge. Further, the 
image of Oq Q Aq is the set of G-parking functions of G. 

Let V = V\{vo}, and let 

X = {veV:{v,vo}^E}. 

Define Ix = Z^^gx ^' ^ configuration on G (having chosen vq as the sink). 
Let K{G)o denote the graph G but with an edge {v,vo} added for each 
V £ X, and fix vq as its sink. Thus, configurations on G and on K{G)q 
are elements of ZF, the free abelian group on V , a subgroup of "LV , the 
configurations on K{G). Recall that K{G) is the graph used to define quasi- 
superstables. 

Theorem 26. 

(1) Define 

(So = {c € (S : c„o = —1 and c + vq> — Ix}- 

Then 

tp{Oo) = {c — vo ■ c + Ix a superstable on K(G)q} 
= So. 

(2) Define 

5o = {c G 5 : c^(j = —1 and c + > 0} C So- 
Then 

tp{Oo) = {c — Vo : c a superstable on G} 

= {c — Vo ■ c + Ix a superstable on K{G)o and c > 0} 

= So. 
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Thus, V(*0o) = ^0 is the set of G-parking functions with respect 
to Vo- 

Proof. We first prove 

Claim A: If c G ZV, then c+lx is superstable on K{G)q if 
and only if c — + ly is superstable on K{G) and c > — Ix- 

For any graph H, let Eh denote its edges. Recall that q is the sink 
vertex for K{G). Let c G IjV with c > — Ixi and consider c — vq + ly as a 
configuration on K{G). Let J7 C y. If G U, then we cannot legally fire U 
from c— uq + Iv: since ic—vo + lv)vo = a-iid {^Oi q} £ Ek(g)j firing J7 would 
result in a configuration with a negative -yo-component. Thus, c — vq + Iv is 
superstable if and only if there are no nonnempty subsets U C.V that can 
be legally fired. 

For V e U QV, consider the edges incident with v that lead out of U: 

M{v, U) = {we{VL} {q}) \U:{v,w}e Ek(g) } 
M{v, U)o = {weV\U -.{v^w} e Ek(g)o } • 

Then, 



M{v, U) \ M{v, U)o = {q} and M{v, U)q \ M{v, U) 
So the cardinality oi M{y,U) is 
(5) \M{v,U)\ = 



iiv^X, 



\M{v,U)o\ + l iiviX, 
\M{v,U)q\ iiveX. 

Now c—?;o+li/ is superstable on if and only if (c— uq+Iv)?; < \M{v,U)\ 

for all w G for all nonempty C y, and c+lx is superstable on K{G)q 
if and only if (c+ lx)v < \M{v, for all -y G ?7 for all nonempty U CV. 
For t; G F we have 

'(c+lx)^, + l iff^X, 



{c-Vo + lv)v = Cy + 1 



{c+lx)v if vex. 



So Claim A follows from (5). The condition c > — Ix is required in the 
statement of the claim since superstables must be nonnegative. 
Now let 

P = {c — Vo ■ c + Ix superstable on K{G)o}. 
The fact that P = So follows directly from Claim A. We now show that 
V'(Oo) = So to finish the proof of part (1). Let O G Oq. Then ip(0) G S 
by Theorem 19. Since vq is a source for O, we have ijj{0)yg = —1, and if 
V G V\X, then ip{0)y > 0. Thus, ijj{0) G So- Conversely, given c G So, run 
the superstables algorithm from Section 3.3 with any vertex ordering of V 
in which vo appears first. Using the notation from the initialization stage 
of the algorithm, let b = Cmax — c. Since c„g = —1, the vertex vq is unstable 
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in b and will fire first. Since c + vq > —Ix, no vertex v & V \ X is unstable 
in b. So when vq fires, the algorithm will orient each edge incident on vq out 
from Vq. If O is the scmiorientation produced by the algorithm, it follows 
that O G Oo and, by Theorem 19, we have ip{0) = c. Thus, ip{Oo) = Sq. 
To prove part (2), let 

N = {c — Vq :ca supcrstable on G}, 

P = {c — Vq : c + Ix a superstable on K{G)q and c > 0} C P. 

Prom part (1), it follows directly that tpi'^o) = P = ^o- To show N = Sq 
and finish, proceed exactly as in the proof of Claim A. Let c G ZV with c > 0. 
We must show that c is superstable on G if and only if c — vq + ly is 
superstable on K{G). Given v dV , this time consider the set 

M{v, U)g = {w^V\U : {v, w} G Eg}, 

and note that \M{v, U)\ = \M{v, U)g\ + 1, with M{v, U) defined as before, 
from which the result follows. □ 

Part 2 of the following corollary recaptures a well-known result from sand- 
pile theory. For instance, it occurs as Lemma 5 in [4] as a statement about 
recurrent configurations, equivalent to our statement in light of Proposi- 
tion 7. 

Corollary 27. Order the superstable configurations on G by the relation 
"<" defined in Section 2.4- Let 5max denote the maximal quasi- super stables 
on G, and let g := — |y| — 1 be the genus of G, as in Corollary 20. Let c 
be a superstable configuration on G. Then, 

(1) c is maximal if and only if c — vq G 5o H iSmax; 

(2) deg(c) := Ylv^v ^ 9 with equality if and only if c is maximal. 

Proof. By Theorem 26 (2), c — G cSq. Hence, part (2) follows immediately 
from part (1) by Corollary 20 (2). We now prove part (1). 
{=^) Suppose c — vq ^ 5niax- By Theorem 26 (2), there exists O G Oq such 
that ip(0) = c — vq. By Corollary 20, we have O ^ Omax, i-e., O is not an 
acyclic orientation of G, but by the beginning of the proof to Corollary 20, 
there exists O' G Omax such that O is a proper subset of O' . It follows that 
V'(O) = c-vq< xp{0') =: c', where c' G S. Since O C O' and O' is acychc, 
we have c^^ = —1 and d + vq > 0, i.e., c' G Sq. Define c' by d = c' — vq. 
Then c' is a superstable configuration on G by Theorem 26, and c < c'. 
So c is not a maximal superstable. 

(<^) Conversely, suppose c — fo G 5max- Take any superstable configu- 
ration c' such that c < c'. Then c' — vq & S by Theorem 26 (2), and 
c — vq < c' — Vq. By maximality of c — vq, it follows that c = c' . Hence, c is 
a maximal superstable. □ 
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Definition 28. The {G,VQ)-semiorder arrangement, denoted J^o, is the set 
of hyperplanes in M" given by 

Xi Xj = 1, 

for all i,j not equal to such that {vi, Vj} G E. 

Definition 29. The regions of J^o, denoted 7?.o, are the connected compo- 
nents of M" \ Jq. 

Define the subset 

Tq = {(xo, . . . , Xn) e M"+^ ■.Xi>XQ + l whenever {vi, vq) G Eq} 

and let 

7^'o = {r G 7^ : r C Tq}. 

The elements of TI'q are exactly those regions with corresponding semiori- 
entations (under p) having vq as a source. Hence, the bijection p: O — > 7^ 
restricts to a bijection Oq — >■ TI'q- The projection mapping (xq, . . . ^ 
(a:;i, . . . omitting the 0-th coordinate, induces a bijection 

TTo : T^-o — > TZq. 

Therefore, we have the following theorem. 

Theorem 30. The mapping 

Po ■= TTo op: Oo -)• Tlo 

is a bijection. 

Example 31. The {G, 'Uo)-semiorder arrangement for graph G of Figure 2 
with Vq = q is drawn in Figure 3. Its regions, TZq, are the projections of the 
regions TZ'q displayed in Figure 6. □ 

The central region of J^o is the region defined by |xj — Xj| < 1 for all 
distinct i,j not equal to such that {vi,Vj} is an edge of G. Inductively 
label the regions of J^o as in Section (3.4), but starting with the central region 
labeled with the configuration that assigns to all Vi such that {vi,VQ} € E 
and —1 to all other vertices, including vq. For each r G TZq, define Ao(r) to 
be the label assigned to r in this fashion. 

Define To = Pq^. There is a version of Theorem 3.4 in this context (proved 
similarly) : 

Theorem 32. We have Xo = ip o tq. So there is a surjective mapping 

Ao : TZq 5^ Sq 

and, for each a € ©O; « commutative diagram 




no, 

where rjafi is the restriction of rj^r to Sq . 
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For i = 1, . . . , n, if {vi,vo} G Eq, let 




otherwise, if {wj, fo} ^ Eq, let 

Q) = {{xi, . . . , Xn) G M" : a;i > + 1 for some j with Vj} G E'}. 

Define = HLi Kf,,o)- 

Definition 33. The admissible regions of J^Oj denoted 7?.o, are the connected 
components of \ ^o- 

Note that TZo = TZo exactly when each nonsink vertex is connected by an 
edge to the sink. 

The admissible regions of are exactly those regions whose correspond- 
ing G-semiorientations satisfy indegQ^Vi) > 1 for all i ^ 0. So combining 
Theorem 26 and Theorem 32 gives 

Theorem 34. 

-^o(^o) = So = {c — vq : c a superstable on G}. 

Example 35. The 9 admissible regions in Figure 3 are labeled by the 8 
distinct superstables (or G-parking functions if one remembers that the sink 
is labeled by —1). The zero-configuration appears twice. □ 

5. Conclusion. 

Let A be an (n+1) x (n+l) matrix. Define Ha to be the set of hyperplanes 

Xi Xj = Aij , 

for all i ^ j such that {vi,Vj} G E. For example, Ha = J' \i A has all Is as 

its entries. 

Define the regions of Ha, denoted IZa, to be the connected components 
of M""*"^ \ Ha- The set of inequalities Xi — Xj < Aij for all i and j defines 
the central region of Ha- We say that Ha has a central region if this central 
region is nonempty. 

Conjecture 36. Suppose that Ha has a central region. Labeling the regions 
of Ha as in Section 3.4 defines a surjection 

Xa '. T^A — ^ S. 

A similar conjecture holds if one first chooses vq as a sink: replace A with 
an n X n matrix, and label regions as in Section 4. So the central region 
would be labeled with the configuration that assigns to vertices connected 
to Vq and —1 to the other vertices (including vq). We conjecture that the 
nonnegative configurations that arise as labels arc exactly the G-parking 
functions. The G-Shi conjecture of Duval, Klivans, and Martin is a special 
case. 
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In the spirit of [1] and [16] , it would be interesting to extend our results to 
the case of multigraphs: graphs in which multiple edges are allowed between 
vertices. If there are k edges between Vi and vj, one might replace the two 
hyperplanes Xj — Xj = ±1 with the 2k hyperplanes Xi — Xj = ±1, . . . , zizk. 
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